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Abstract 

In the study of the small ten-dimensional Schwarzschild blackhole, the blackhole to 
string transition is an important problem. In [1], a possible identification is made 
between the Gross- Witten-Wadia (GWW) type third-order large- N phase transition 
in the boundary gauge theory and the string-black hole transition in the bulk. In this 
paper, we exhibit the existence of the GWW transition by Monte Carlo simulation 
in the zero mode bosonic action of the finite-temperature Af = 4 SYM theory on 
S 3 . Exhibiting this transition in the truncated but highly non-trivial gauge theory 
implies that in the vicinity of the critical temperature T c , the system goes critical, 
and the fluctuations give rise to universal formulas derived in [1]. We also discuss 
the issue of SO (6) R-symmetry breaking. 
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1 Introduction 

Understanding the string-black hole transition is an important problem in string theory. 
The radius of a Schwarzschild blackhole becomes smaller with rising temperature and at 
a certain temperature the blackhole transits to a gas of strings [2-6]. This is a difficult 
problem to address as it needs a proper understanding of non-perturbative effects in string 
theory. The finite temperature of a Schwarzschild blackhole breaks supersymmetry and 
string loop corrections are uncontrolled. In [1] it was demonstrated that the problem 
of the string-blackhole transition can be formulated in a space-time with AdS^ x S 5 
boundary conditions. This enabled one to use the AdS/CFT correspondence [7] to map 
the string-blackhole transition phenomenon to a Gross- Witten-Wadia (GWW) type phase 
transition [8-10] in the boundary gauge theory defined on S 3 x R. Q 

The compactness and positive curvature of the space S 3 permits one to integrate out 
all other modes to get an effective multi-trace unitary matrix model for the zero mode 
of the Polyakov line. Based on the works [11-18], this type of effective unitary matrix 
model was analyzed in [1] to show the existence of the GWW type transition. The o(l) 
part of the gauge theory effective action was also calculated in a double scaled region near 
the transition temperature. The o(l) part is universally given in terms of F(t), where F 
satisfies the following differential equation 



and f(t) is the Painleve II function, and t is a scaled variable proportional to (T — T C )N3. 

The derivation of the effective unitary matrix model from the gauge theory is a subtle 
one. In a weakly coupled gauge theory one may demonstrate this explicitly in perturbation 

4 One should caution against too literal an interpretation of the gauge theory result since at the cross 
over, the winding Polyakov line is non-zero, signaling that a black hole (without a space-time description, 
still persists at the phase transition and passes over into a stringy description only beyond the temperature 
at which the phase transition occurs. We would like to thank Juan Maldacena for a discussion of this 
point. 
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theory at large N [19]. However the situation is less clear in the strong coupling regime. 
One difficulty comes from the Gregory-Laflamme transition for a small AdS$ x S 5 black- 
hole. This transition breaks the SO (6) symmetry of S 5 and the question arises whether 
there is a new zero mode associated with this transition, and whether the unitary matrix 
model is a good description after this transition. In [1] it has been shown, using a super- 
gravity analysis within the AdS/CFT correspondence, that even at strong coupling, the 
unitary matrix model serves as an effective description. 

Given the physical relevance of the GWW transition, it is important to see if this 
phenomenon occurs when one is not dealing with the dynamics (however complicated) 
of a single unitary matrix or the quantum mechanics of a single unitary matrix. It is 
not at all obvious that this large- N transition occurs in more complicated models of non- 
commuting matrices and gauge theories. In the past this question has been explored by 
Douglas and Kazakov [20] in their study of two-dimensional Yang-Mills theory on S 2 . 
However this problem too, gets recast into a problem of a single unitary matrix because 
the partition function turns out to be the heat kernel on the unitary group. 

In order to answer these questions, there seems to be no analytic tools as is usually 
the case with complicated dynamical problems. Hence we use numerical Monte Carlo 
methods to explore and exhibit the large- N GWW transition and also study the question 
of R-symmetry breaking at large N. Since the full M = 4 SYM theory on S 3 x R is 
too difficult, in the first run we study the gauge theory restricted to the zero modes of 
the bosonic sector. It is likely that this reduction captures the essential features of the 
dynamics. It is motivated by the fact that the metric of the small Schwarzschild black 
hole is uniform on S 3 . Regarding fluctuations in the bulk, the zero mode gauge theory 
has correspondence with fluctuations in the bulk which are independent of S 3 and only 
depend on the radial AdS§ coordinate and time. 

The importance of exhibiting this transition lies in the fact that in the vicinity of the 
critical temperature T c , the system goes critical and the fluctuations give rise to universal 
formulas (II. ip which solely depend on the multi-critical point which is characterized by 
the exponent | in the scaling law (T — T c ) ~ N~» . Hence the formulas for the black hole 
cross over which were derived using the effective unitary matrix model in [1] are also valid 
while working directly with the zero mode sector of the gauge theory. 

This paper is organized as follows. In Section 2, we introduce the zero mode action 
of the bosonic part of the M = 4 SYM theory on S 3 x R. In Section 3, we discuss the 
numerical studies of the GWW-type phase transition. In Section 4, we study the SO (6) 
R-symmetry using Monte Carlo simulation. Section 5 is devoted to conclusions and future 
outlook. 



2 The model: J\f = 4 SYM theory reduced on S 3 



We study the M = 4 SYM theory, when all the bosonic fields are restricted to their zero 
modes on S 3 . 




(2.1) 
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(2.2) 



and D t is a covariant derivative defined by 

DtM^t) = dtM^t) - i[A(t), M^t)]. (2.3) 

D is the dimensionality of the model, and the dynamical variables A(t) and M^(t) (fi = 
1, 2, • • ■ , D) are N x N Hermitian matrices, which can be regarded as the gauge field and 
the SO(D) adjoint scalars, respectively. This model has a U(N) gauge symmetry 

M,(t) -> g(t)M,(t)g\t), A(t) -> ^(t)^) + ig{t)^-. (2.4) 

The Euclidean time t in the action f)2.2p has a finite extent (3, which is the inverse tem- 
perature (3 = 1/T. Both the gauge and the scalar fields obey the periodic boundary 
conditions 

A(t + /3) = A(t), M lt (t + p) = M lt (t). (2.5) 

While this model has three parameters, j3, A and m, these are not independent of each 
other, as m can always be set to unity by the following redefinitions 

/*--, (2.6) 
m m 6 

and rescaling of the fields 

A(t)^-A(t), M^t)^m^M^t). (2.7) 
m 

The periodic boundary condition (12.51) prevents us from fixing the A = gauge. However 
we can fix a gauge where the gauge field is static and diagonal: 

A= -diag(a 1 ,a 2 ,---,ajv), (2.8) 

where a a G (— 7r,7r]. The indices a, 6, ■ • • run over 1,2, ■ • • , iV. This gives rise to the 
Fadeev-Popov term 

N 

St P . = - logsin|K-a 6 )/2|, (2.9) 

a,b=l,a^b 

whose derivation is given in full detail in [21,22]. In the following, we study the action 

S = S' + S Lp .. (2.10) 

We study the model (12.10 1 ) numerically by Monte Carlo simulation. The details of the 
algorithm are given in [23]j. We simulate the model with the time direction discretized. 
We apply the heat bath algorithm to the scalar fields and Metropolis algorithm to the 
gauge field, respectively. It turns out that taking 10 lattice points of the time direction is 
enough and that increasing the lattice points further does not affect the result. 



5 In [24,25], they discuss the simulation of the supersymmetric gauge theory at finite temperature 
using Fourier transformation, instead of lattice discretization. 
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3 GWW phase transition 



In this section, we study the GWW type phase transition of the simplest unitary matrix 
model, for which an analytical solution is available [8-10]. This model is useful to test 
the accuracy of the numerical method. 



3.1 The trW + trU model 

We start with the numerical simulation of the unitary matrix model consisting only of 
trU without adjoint scalar fields, where 

U = V exp [ % f fi dtA(t)) . (3.1) 



V denotes the path-ordered product. We consider the partition function 

Z g = J dU exp (^{trU + tr Eft) J . (3.2) 

and define u n = jjtrU n for an integer n. In the static and diagonal gauge (12. 8p . 

1 N 



N I 



These are the moments of the density of eigenvalues: u(a) = -k- Yla=i 
The first two moments are given by 



(Ki) 
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(9 < 1) 



The third-order transition at the point g — 1 is the GWW transition. This is a 
transition between the gapped and ungapped phases of the eigenvalue distribution of the 
unitary matrix model. For a generic unitary matrix model, all u n 's show a similar non- 
analytical behavior like ui, because near the gap opening point, the relevant operator is 
given by a linear combination of u n [1]. 

We first verify this result numerically using Monte Carlo simulation. To this end, we 
take static and diagonal gauge (12. 8p and add the Fadeev- Popov term (12. 9p . Namely, we 
apply the Metropolis algorithm to the action 

_^( tr E/ + trE/t)- logsin|^p^|. (3.5) 

We plot the VEV's (|wi,2|) against g in figure [TJfor N = 128, and find that they actually 
agree with the result (13. 4p . 
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Figure 1: The vacuum expectation values (|ni i2 |) against g for iV = 128. 

3.2 GWW phase transition in the gauge theory reduced on 

We next study the saddle point of the gauge field by adding the chemical potential /i(tr U+ 
trW) to the action (12.101) . Namely, we study the matrix model 



where the terms S' and 5 g .f. are defined in ( 12. 2ft and ( 12.91) . respectively and U is the 
Polyakov line defined in ( 13. ip . 

3.2.1 D = 2 case 

We first study the D = 2 case, in which the numerical simulation of large N is reachable 
at a reasonable CPU time. The phase transition of the one- dimensional matrix quantum 
mechanics with respect to the temperature has been studied in [23,26-28] in the absence 
of the chemical potential. The Polyakov line (\ui\) is small in the low-temperature region, 
while it is large in the high-temperature region. We focus on the low-temperature region 
(3 = 2.0, in which the Polyakov line (|iti|) is small for il = 0.0. We plot the result of the 
D — 2, A = m — 1.0 and N = 48 case in figure [2j 

The graph above indicates a signature of the phase transition (possibly third or higher 
order) near the critical point 



at which the Polyakov line is = 0.5. This is expected from the fact that (|ii2|) ^ 

(| «i |) near the transition point. To understand the nature of the transition we first numer- 
ically plot the derivative d ^ 1 " in figure [2] (left). The derivative seems to be continuous, 
and hence the possible transition should at least be of third order. Numerical errors pre- 
vent us from going further and calculating the higher derivatives directly from our data. 
Instead, in figure [2] (right) we try to fit our data with analytic functions in the regime 
H < fi c and fi > fx c and extrapolate the information about derivatives from the fitted func- 
tions. It should be noted that the fitted functions do not necessarily represent the correct 
analytic form of the exact answer, but they can be viewed as a close approximation. 



S g = S' + 5 g . f . + N0n(tT U + tr U ] ), 



(3.6) 



He ^ 0.22, 



(3.7) 
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We fit the VEV with the function 



We exploit the fact that in the large-iV limit, (|iti|) is at /i = and -> 1 as 

dfi 



ji — > oo. And from the fact that and its first derivative ^j" 1 ^ are continuous at the 



critical point /i = /i c , we obtain the condition for r\ and r 2 . 

The parameters qi^ are fitted as 

gi = 0.503542 ±0.01181, q 2 = 0.53791 ± 0.003644. (3.10) 

In this case, the coefficients are T\ = —0.0121 and r 2 = —0.029318, respectively. We 
find that the contribution of the terms ^i(^j-) 2 for fi < fi c and r 2 (^j-)~ 2 for fi > \x c is 

small compared to the rest of the terms in (13. 8p . Since its second derivative d is 
discontinuous at the critical point /i = /i c , this system undergoes the GWW type third- 
order phase transition. 

The VEV ( 1 ) is small when fi < fi c , and in this region (|w 2 |) is closer to zero at 
larger N. When fi > fi c , it is fitted with the following function similarly to the unitary 
matrix model. 

(M> = i-^ + 4. (a*>aO- (3-n) 

fi fi z 

Next, we turn our attention to the high-temperature case f3 = 0.2, in which the VEV's 
of the Polyakov line ( |^i,2 1 ) are large even in the absence of the chemical potential. We 
plot (|w 1)2 |) against \i for D = 2, (3 — 0.2, A = m = 1.0 case for iV = 48 in figure 
[3j In contrast to the low-temperature case, we find that there is no GWW type third- 
order phase transition in this case and that the VEV's of the Polyakov line increase 
monotonically. 
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Figure 3: The vacuum expectation values ( |^i,2 1 ) against fi in the D = 2, (3 = 0.2, 
A = m — 1.0 case for N = 48. 

3.2.2 D = 6 case 

We next study a different dimensionality, D = 6. Similarly, we plot the VEV's { |^i,2 1 ) 
against /i in the D = 6, (3 = 2.0, A = m = 1.0 case for = 16 in figure HI We read off 
the critical point as 

fi c ~ 0.20. (3.12) 

Then, we fit them with the functions (13. 8p and (13.111) . In this case, the parameters are 

q x = 0.53803 ± 0.02015, q 2 = 0.542291 ± 0.006888. (3.13) 

The coefficients r 12 are r\ = —0.039095 and r 2 = —0.0412257, which suggests that the 
contribution of the r 12 terms is smaller than that of the rest of the terms in (13.81) . We 
find that the result is similar to the D = 2 case. 
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Figure 4: The vacuum expectation values ( |^i,2 1 ) against fi in the D = 6, (3 = 2.0, 
A = m = 1.0 case for A^ = 16. 
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4 SO (6) R-symmetry breaking 



In this section, we study the spontaneous breaking of the SO (6) R-symmetry of the 
model f)2.10p by Monte Carlo simulation. Throughout this section, we focus on the D = 6 
dimensional case. In analogy with the IKKT-type matrix model [29], we consider the 
following observable [30-36], 

r iv (t) = jjti{M lt {t)M v {t)). (4.1) 

In our case, we integrate the operator I'„ v {t) with respect to the time direction and obtain 
the "integrated moment of inertia tensor" 

"P 



1 r 

Iflu = Nj dUrM ^ t ) M ^ t )- ( 4 - 2 ) 



We define the eigenvalues of this 6x6 matrix 1^, which are all real positive, as A M with 
the specific order 

Ai > A 2 > ■ ■ • > A 6 . (4.3) 
We consider the following 50 (6) invariant quantity [30], 

J = - ( 4 - 4 ) 

This quantity measures the variance of the eigenvalue distribution of I^ v . Using large- N 
factorization we get 

(J) = (^W-(^) 2 ) (4-5) 



(4.6) 



Using the fact that the VEV of any 5*0(6) two-tensor is proportional to 5^, i.e. (I^ v ) = 
A^jy, we get (J) = 0. This relationship is not true in general and we expect (J) ^ 
at finite N. In the case when (J) is non-zero, the width of the eigenvalue distribution 
of 1^ is non-zero. The above scenario implies that the dominant contribution of the 
path integral comes from the configurations for which the eigenvalues of are not equal 
and consequently the SO(6) symmetry is broken. Hence by plotting the VEV's of the 
eigenvalues of and measuring the width of the distribution, we can figure out the 
possibility of 5*0(6) symmetry breaking at large N . This leads us to evaluate the VEV's 
of these eigenvalues (A M ) in the large- N limit. After diagonalization, the residual 50(6) 
transformations permute the eigenvalues A M . Hence an unbroken 50(6) symmetry implies, 

(A 1 ) = (A 2 ) = --- = (A 6 ). (4.7) 

whereas a broken 50(6) symmetry implies that for some \i > u, (A^) > (\ u )- To this end, 
we extrapolate the large- N limit from the simulation of finite N. If the eigenvalues (A^) 
are all equal in the large- N limit, this suggests that the 50(6) symmetry is unbroken. 
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4.1 Dynamical gauge field 



We first study the SO (6) R-symmetry breaking when gauge field A is integrated. To this 
end, we update the gauge field A, as well as the scalar fields M^{t) via the usual algorithm. 
We extrapolate the large- N limit from the finite- N results of N — 16, 20, 24, 28, 32 for 
the high-temperature (3 = 0.1 and the middle-temperature (3 = 1.0 cases. We plot the 
eigenvalues (A M ) against 1/N in figure [5) 
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1/N 




0.01 0.02 0.03 0.04 0.05 0.06 0.07 
1/N 




0.01 0.02 0.03 0.04 0.05 0.06 0.07 
1/N 



0.01 0.02 0.03 0.04 0.05 0.06 0.07 
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Figure 5: (A M ) against 1/N for m = 1.0, and N = 16,20,24,28,32. (/3, A) = (0.1,0.1) 
(upper left), (/?, A) = (1.0,0.1) (upper right), (/?, A) = (0.1,1.0) (lower left), (/3, A) = 
(1.0, 1.0) (lower right). 

It turns out that the eigenvalues (A M ) converge to the same value in the large- N limit. 
This behavior is qualitatively the same for other parameter regions of the action ( 12.101) . 
This indicates that the 50(6) R-symmetry is unbroken in the matrix model (I2.10p . 

4.2 Uniform and Clumped configurations of the gauge field 

We next study the 50(6) R-symmetry breaking in the specific configurations of the gauge 
field, which correspond to AdS 5 x S 5 and a black hole. To this end, we put the constraints 
on the gauge fields. In the following, we focus on the high-temperature j3 = 0.1 and 
massive m = 1.0 case, and take A = 1.0. 
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1. Uniform distribution 



We take the diagonal part of the gauge field (12. 8p as 



7T 



(2a- N). 



(4.8) 




In this case, the Polyakov line U satisfies (\u n \) = for any nonzero integer n. In 
the AdS/CFT correspondence uniform distribution, which is depicted in figure E] (1), 
corresponds to the AdS 5 x S 5 spacetime [1,37]. To realize this configuration, we skip the 
Metropolis algorithm to update the gauge field A and fix the configuration of the gauge 
field to be (14.81) . We update only the scalar fields M^it) via heat bath algorithm. 

2. Clumped distribution 

In the clumped distribution, we constrain the gauge fields in a small region a a G 
[— 7re, +vre], which is opposite to the gapped distribution. This distribution is depicted in 
figure (2). Similarly to the gapped distribution, we take e = 0.05. If a a goes out of the 
region [— 7re, +7re], we automatically reject that configuration. This configuration coming 
from a gapped distribution of eigenvalues corresponds to the blackhole state as can be 
indicated by an analysis of large- N perturbation theory around the gapped phase [1]. 



Figure 6: Distribution of the diagonal part of the gauge fields {e t/3aa } in (l)uniform and 
(2) clumped distribution. 

Similar to the case when we updated the configuration of the gauge field, we make a 
large- N extrapolation of the eigenvalues (A^ t ). 

We plot in figure [7]the eigenvalues (A M ) against 1/N for the high-temperature (3 = 0.1 
and A = m = 1.0 case. In these cases, too, the eigenvalues (A M ) converge to the same 
value at large N. We find that the 5*0(6) R-symmetry of the scalar field is unbroken for 
these configurations of the gauge field. 

5 Conclusions 

In this paper, we have exhibited the GWW large- N phase transition using Monte Carlo 
simulation in the zero mode reduction of the bosonic part of the TV = 4 SYM theory on 
the S 3 x R space. We have studied the saddle point by adding a chemical potential to 
the reduced action, and observed a third-order phase transition in the large- N limit. Its 
significance is that the large- N transition signals critical behavior and the the properties of 




(l)uniform 



(2)clumped 
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Figure 7: (A^) against 1/N for /3 = 0.1 and A = m — 1.0 in the uniform (left) and clumped 
(right) distribution. 

the model in the vicinity of the critical point are universal. Hence we expect that the o(l) 
free energy is given by (11.11) . We have also numerically found that the 50(6) R-symmetry 
is NOT spontaneously broken, in the large- N limit. In the d = and d = 1 unitary matrix 
models the physical mechanism for the GWW transition is well understood. In the d = 
models the repulsion between eigenvalues, from the measure, and their attraction in the 
potential well, are competing effects which lead to this transition [8,9]J§ In the d — 1 
models the phase transition is signaled when the Fermi level reaches the hump (maximum) 
of the potential [10]. In the more complicated models we have explored, there are typically 
non-commuting matrices and explanation seems to be less obvious. 

In the future it would be instructive to go beyond the zero mode approximation and 
develop numerical methods to include the variation of the fields on S 3 . Also, it would be 
interesting to be able to include the adjoint fermions of the gauge theory in the numerical 
calculation. 
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